Novel glassy behavior in a ferromagnetic p-spin model 
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Recent work has suggested the existence of glassy behavior in a ferromagnetic model with a four- 
spin interaction. Motivated by these findings, we have studied the dynamics of this model using 
Monte Carlo simulations with particular attention being paid to two-time quantities. We find that 
the system shares many features in common with glass forming liquids. In particular, the model 
exhibits: (i) a very long-lived metastable state, (ii) autocorrelation functions that show stretched 
exponential relaxation, (iii) a non-equilibrium timescale that appears to diverge at a well defined 
temperature, and (iv) low temperature aging behaviour characteristic of glasses. 



I. INTRODUCTION 

Many liquids when cooled below the melting temper- 
ature do not crystallize. Instead, as the temperature is 
lowered, the relaxation time (or the viscosity) increases 
dramatically, eventually becoming so large that the liquid 
appears frozen on experimental timescales. For all practi- 
cal purposes the system has solidified yet there is no long 
range order. The system is said to have become a glass. 
Despite over fifty years of work, a consensus is lacking on 
a theory of this phenomenon jl| . Even as basic a ques- 
tion as the existence of a genuine thermodynamic glass 
transition is difficult to decide experimentally. Therefore, 
much work in this field has been numerical. Simulations 
in the continuum of glassy systems remain difficult, and, 
despite the substantial progress in this direction Q , there 
are still many unresolved issues. 

Because of these problems it would be useful to have 
simplified lattice models that behave in some respects 
like glasses and which do not have quenched disorder. 
Recently, Lipowski || and Lipowski and Johnston [|]|| 
studied such a model by Monte Carlo simulations and 
provided evidence of glassy behavior. The model is de- 
fined by the Hamiltonian 

H = -jY,S i S j S k S h (I) 

□ 

where the spins sit on a cubic lattice and the interactions 
are between the four spins at the corners of each plaque- 
tte. The spin variables take values Si = ±1 and J is the 
strength of the ferromagnetic coupling. (We will take 
J = 1 in all that follows.) Their study was motivated 
by the fact that the spin-glass version of this model was 
found numerically || to have several features in common 
with structural glasses. In addition, there is a close con- 
nection H between the plaquette model, Eq. [l], and the 



Ising model with competing interactions studied in the 
context of glasses by Shore and Sethna . Lipowski and 
Johnston ||] showed that, in equilibrium, the plaquette 
model has a first order transition as the temperature is 
lowered but that this transition is never observed in sim- 
ulations within accessible timescales. Instead, one sees 
what appears to be a dynamical transition to a glassy 
state 0,§. 

The work of these authors concentrated on single time 
quantities. To see clear signatures of glassy behaviour 
and to understand which particular aspects of glasses are 
present in this model, it is important to study two time 
quantities such as autocorrelation and response func- 
tions. We do just that in this paper and our findings 
are organized as follows: in section II, we summarize 
the results of Lipowski and Johnston on the statics of 
the plaquette model and the strong metastability that 
they observed. In section III, we examine the metastable 
state in greater detail. In particular, we calculate the 
timescale for nucleation and show that it is indeed ex- 
tremely large. In section IV, we present results for the 
energy auto-correlation function in the supercooled phase 
which is unambiguously described by a stretched expo- 
nential. We investigate the behavior of the associated 
timescale as a function of temperature, showing that it 
diverges with a power law at a temperature indistinguish- 
able from the lower limit of metastability. This temper- 
ature is that which the authors of [||,[| called the glass 
transition temperature T g . In section V, we study the 
dynamical behavior below T g . The spin-spin autocorre- 
lation function and overlap distribution function exhibit 
aging of a form that is characteristic of glasses, show- 
ing that the plaquette model has a complex free-energy 
landscape. These findings provide further evidence to 
support the conjecture that the plaquette model exhibits 
an ideal glass transition. 
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II. THE PLAQUETTE MODEL 

In equilibrium, the model defined by Eq. |l| has a first 
order transition at T c ~ 3.6 [|- This system pos- 
sesses a novel symmetry: nipping all the spins in any 
plane of the cubic lattice leaves the Hamiltonian invari- 
ant, leading to a degeneracy of ground states ~ 2 dL for 
a cubic system of size L in d dimensions. Note, however, 
that the ground state entropy is finite. Because of this 
symmetry, the magnetisation is not a good order param- 
eter. Instead, one can use the internal energy per spin, 
E, to investigate this model further. 

Lipowski and Johnston looked at the dynamics of E 
following a quench from a high temperature state to 
a temperature T < T c [Q. They found different be- 
haviour depending on the value of T. In particular, 
for 3.4 < T < T c , E appears not to relax to its value 
in the low temperature phase. Rather it settles into a 
plateau at an energy higher than that of the true low 
temperature phase at the same temperature. Further- 
more, the time the system spends on the plateau (in a 
"supercooled" state) appears to increase at least expo- 
nentially with system size, suggesting the existence of a 
genuine metastable phase. Such behaviour , if present, 
would indeed be strange because one does not expect 
true metastability (i.e. two local minima of the free en- 
ergy separated by diverging free energy barriers) in sys- 
tems with finite range interactions. Moreover, nucleation 
theory tells us that the larger the system, the shorter the 
lifetime of a metastable state, provided that the system is 
large enough to support a nucleating droplet. To under- 
stand whether the strong metastability is indeed a finite 
size artifact, Lipowski and Johnston studied a quench 
starting from a configuration with a droplet of the low 
temperature phase in a high temperature background. 
They reported that for sufficiently large seeds, the system 
always relaxed to the low temperature phase, concluding 
that ordinary nucleation was operative. However, this 
does not explain why the supercooled phase appears so 
stable on observable timescales. The analysis we report 
in the next section provides an explanation for this phe- 
nomenon. 

For quenches to below T — 3.4 the behaviour was 
found to be quite different. The plateau disappeared and, 
after initial transients, the energy decreased extremely 
slowly. Moreover, it appeared to relax not to the value 
in the low temperature phase, but to a value extensively 
higher. They identified this state as the glassy phase. 
Finally, they have recently shown that a characteristic 
zero-temperature length increases very slowly with the 
inverse cooling rate, as is seen in ordinary glasses j5|. 

To summarize, this model has many properties that are 
similar to those seen in real glasses: a first order tran- 
sition analogous to crystal melting, strong metastability 
effects in the supercooled phase, and a kinetic transition 
to a glassy state characterised by slow dynamics. 

In the next section we will discuss the anomalous 



metastability described above before looking in some de- 
tail at the dynamical behavior of this model by study- 
ing the autocorrelation functions both in the supercooled 
phase and in the "glassy" state. 



III. METASTABILITY 

We begin by summarizing the basics of nucleation the- 
ory JlO| . The idea is that if the system is in a metastable 
state, the process of finding the true free-energy mini- 
mum involves the formation of a droplet of the stable 
phase which can grow and take over the whole system. 
Phenomenologically, one can express the free energy bar- 
rier to forming a droplet of linear dimension R as A = 
AaR d ~ l — BSfR d in d dimensions. Here a is the surface 
tension between the metastable and stable phases and Sf 
is the difference in the bulk free energy densities of the 
two phases. A and B are constants characterizing the ge- 
ometrical shape of the droplet. Because the surface term 
is a cost and the bulk term is a gain, one can maximize 
A to find a critical radius R* = (A(d - \) / Bd){a / 5f). 
Droplets of a radius greater than R* can grow and take 
over the whole systems while those of radius less that 
R* will shrink. At low temperatures and in three dimen- 
sions, the timescale for nucleation of a droplet is then 
given by r nuc - exp(4^ 3 CT 3 / '27 B 2 Sf 2 T). 

For the plaquette problem, the free energy difference 
between the two phases can be estimated from the data 
in 0]. We find Sf « 0.5 \T - T c \. The surface tension 
a at T c can be calculated from the size dependence of 
the mean time it takes the system to flip between the 
high and low temperature phases. This time Tfn p is given 

by Tflip ~ e 2aL / Tc , where L is the system size and the 
factor of 2 comes from using period boundary conditions. 
We have performed standard single spin-flip Monte Carlo 
simulations to determine this transition time numerically 
for systems of size L = 4, 5, 6, 7, 8. The transitions are 
identified from jumps in the energy between the values 
appropriate to the high and low temperature phases. For 
each system size, an 'effective transition temperature', 
T C (L), is determined by the requirement that the system 
spend equal time in both phases, i.e. at each of the two 
possible values of the energy. This gives T C (L) ~ 4.45, 
4.19, 4.03, 3.92 and 3.86 for L = 4, 5, 6, 7 and 8 respec- 
tively. The mean transition time, Tfl; p (L), is measured 
at T C (L), and the surface tension, a, obtained by plot- 
ting lnTflip(L) against L 2 /T C (L). The data is presented 
in Fig. |l|. For large L, the points should approach a 
straight line with slope 2a. Since there is clear curvature 
in the data for the small values of L available, we attempt 
an extrapolation to large L by plotting [T C (L)/L 2 ] Inrgip 
against T C L/L 2 , as shown in the inset to Fig. The in- 
tercept on the vertical axis gives 2a. While the extrapo- 
lation to L — oo is necessarily subjective, from the curva- 
ture of the data it is clear that any reasonable procedure 
will give a value of 2a greater than the value ~ 0.62 ob- 
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tained from a linear extrapolation of the last two points, 
and less than the value ~ 0.69 given by the L — 8 point. 
To obtain lower bounds on critical drop sizes and nucle- 
ation times we use the estimate a ~ 0.31 obtained from 
the linear extrapolation. Assuming that a depends only 
weakly on T close to T Cl and taking T c — 3.6, we estimate 
the critical droplet radius at T = 3.5 to be R* w 25 for 
a cubic droplet (A = 6, B = 1) and rj 12 for a spherical 
droplet (A = An, B = Att/3). The corresponding time to 
nucleate a cubic droplet is of the order of 10 47 MCS (and 
- 10 25 MCS for a spherical droplet). 
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FIG. 1. Logarithm of the flip time between high and low 
temperature phases at T c as a function of system size. The 
errors are smaller than the symbols plotted. The inset shows 
the extrapolation for large L. 

These findings provide an explanation for the results 
of Q], in which it was noted that the relaxation time of 
the metastable phase grows rapidly with system size. In 
||, the system sizes simulated were all less than, or of 
the same order as, R*. In this limit, the free energy bar- 
rier to nucleation is set by the system size alone, as a 
droplet of the equilibrium phase cannot fully form. Con- 
sequentlythe lifetime of the metastable state will grow 
exponentially with the size of the system [0] as reported. 
We thus conclude that the supercooled state is indeed 
metastable but its lifetime is many orders of magnitude 
larger than the times accessible in Monte Carlo simula- 
tions. 



IV. DYNAMICS IN THE SUPERCOOLED PHASE 



The philosophy underlying the search for models with- 
out disorder that show glassy behavior is that in glasses, 
inhomogeneities are self-generated, even though the mi- 
croscopic Hamiltonian has no intrinsic disorder. Thus 
there is a spontaneous separation of the degrees of free- 



dom into variables that exhibit either slow or fast relax- 
ation. One expects that it is the slow degrees of freedom 
that freeze as the glass transition is approached. How- 
ever, this separation is not necessarily a priori obvious, 
and one has to be careful about which quantities one 
studies. It might happen, for example, that the spin-spin 
autocorrelations do not show any anomalous behavior, 
but that some complicated functions of the spins show a 
dramatic slowing down. For this reason it is often useful 
to look at autocorrelation functions of global quantities 
because these presumably pick out the longest timescale 
or the slowest degree of freedom. 

For the system defined by Eq. |Tj, one such quantity is 
the energy per spin, E, and we now discuss our findings 
for the energy autocorrelation function in the range of 
temperatures T g < T < T c . Our results have been ob- 
tained from Monte Carlo simulations using single spin- 
flip Glauber dynamics. We have considered cubic sys- 
tems of linear size L, with L ranging from 12 to 64. 
Most of the data presented here is for a system size of 
L = 48. To measure the energy autocorrelation function 
at some temperature T we start from a random configu- 
ration (infinite temperature) and quench to the temper- 
ature T . We then wait for a time t w and measure the 
autocorrelation function A(t,t w ) for subsequent times, 
t, where A(t,t w ) =< E(t w )E(t w + t) >. The angular 
brackets denote averages over random initial conditions. 
In equilibrium one expects A(t,t w ) to be independent 
of t w and it is only in equilibrium that one can define 
a meaningful timescale associated with relaxation. In 
the supercooled phase, we found that equilibrium is eas- 
ily achieved in terms of this definition even though the 
system is in a metastable state. In particular, for any 
temperature T g < T < T c , the autocorrelation function 
becomes independent of t w once t w becomes of the order 
of 10 3 MCS. Note that this time is at least 10 orders of 
magnitude smaller than the lifetime of this metastable 
state. 

In Fig. H we show A{t,t w ) for three different temper- 
atures below T c and above T s , together with attempted 
fits to stretched exponentials of the form 

A(t,t w )=Aoe-^y. (2) 

As can be seen from the figure, the fits are extremely 
good and we can identify a timescale for relaxation, r, at 
each temperature. The exponent [3 is also temperature 
dependent and is found to lie in the range 0.5 < p < 1. 
In Fig. U we plot the behaviour of r as a function of T. 
It appears to increase sharply as T is reduced and can be 
fit quite accurately by a power-law divergence , 

r = ^^, (3) 
T-3.39 v ' 

as shown in the figure. This gives a "critical" tempera- 
ture very close to that which was identified in ||] as the 
lower limit of stability, T g , and interpreted as the glass 
transition temperature. Both stretched exponential re- 
laxation functions and a diverging relaxation time are 
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features of real glasses. We also tried a Vogel-Fulcher 
(stretched-exponential) fit to r but found that the cor- 
responding "critical temperature" is significantly lower 
than T g . As we will show in the next section, T g signifies 
the onset of aging phenomena and it is thus impossible 
to define a timescale below this temperature. (We have 
also measured the spin-spin autocorrelation function but 
found that the behaviour was much harder to interpret. 
In fact, the difficulty of interpreting relaxational behav- 
ior of local quantities has been noted before in molecular 
dynamics studies of Lennard- Jones systems in two di- 
mensions 
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FIG. 2. The energy autocorrelation function A(t,t w ) in 
the supercooled phase. Three temperatures are shown: 
T = 3.6(o), 3.5(A), 3.42(o). The lines are fits to stretched 
exponentials, Eq. H, with parameters: t = 11,(3 — 0.7 at 
T = 3.6, t = 21,(3 = 0.66 at T = 3.5 and t = 12,(3 = 0.6 at 
T = 3.42. 

A possible explanation for this divergence is to note 
that in mean-field theory, the plaquette model has a first 
order transition at a temperature quite close to T c Q . In 
any mean-field first order transition, one expects a region 
of metastability delimited by the spinodal temperature. 
Above this temperature, there is genuine metastability 
and one expects timescales of autocorrelation functions 
to diverge as the spinodal temperature is approached 
When one moves away from mean-field, genuine metasta- 
bility is no longer present and instead the spinodal gets 
rounded out. Our results for the autocorrelation func- 
tions, and those of the authors of m for the energy re- 
laxation, suggest that the spinodal is actually remarkably 
well-defined in this short-range system. It is therefore 
plausible that what we are calling the glass transition 
temperature is really a remnant of the mean-field spin- 
odal. This too is something that is seen in Lennard-Jones 
systems in two dimensions pT| . However, the stretched 
exponential form of A(t, t w jwould not be expected in 
any simple mean- field treatment. 

Thus, the short-range plaquette model has a first or- 
der transition analogous to crystal melting, but because 



of the extraordinarily long-lived supercooled phase, it 
also exhibits timescales for relaxation that diverge at 
T g . We would like to emphasize that in this model T g 
is not strictly well-defined. In fact, we expect that as one 
gets closer and closer to T g , or waits for longer in the 
metastable state, the divergence we are seeing will get 
rounded out. Numerically, however, it is not possible to 
see this rounding in any reasonable amount of time. This 
again is a feature in common with glass forming liquids. 
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FIG. 3. The timescale of relaxation r as a function of tem- 
perature. The curve is the best fit to a power-law divergence 
given by Eq. H. 



V. DYNAMICS IN THE GLASSY PHASE 

For a quench to a temperature T > T g , as we have seen 
above, the system becomes trapped in a metastable state 
for times up to I0 24 MCS. Within this time window, two- 
time correlation functions become waiting-time indepen- 
dent and the system behaves as if in equilibrium. This 
is no longer the case if T < T g . The system now contin- 
ues to evolve slowly for all accessible times and two-time 
quantities exhibit aging. To investigate this behaviour we 
have measured the spin-spin autocorrelation function, 



C(i, t w ) — — ^ ^ Si(t w ) Si(t + t w ), 



(4) 



for a range of waiting-times t w . Our results are shown in 
Fig. H There is a strong dependence on t w for all times 
and temperatures we have considered. In the long time 
regime, C'(t, t w ) decays as a power of t with an exponent 
of around 0.35. Similar behaviour has been noted in MD 
simulations of binary Lennard-Jones mixtures Jl^]. We 
have also attempted to collapse the data at a fixed tem- 
perature with a scaling form C(t, t w ) — C(t/T(t w )) where 
r(t w ) ~ t^. While the fit is not perfect, our best estimate 
for a is a — 2. This type of behaviour has been termed 
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super-aging |51| as the characteristic timescale for relax- 
ation, r(t w ), grows faster than the weighting time itself. 
Note that the spin-glass version of this model exhibits 
sub- aging with a w 0.77 0. 
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FIG. 4. The spin autocorrelation function C(t,t w ) in the 
glassy phase at temperature T — 3.2. Eight different different 
waiting times are shown. From bottom to top t w — 10, 100, 
200, 500, 1000, 2000, 5000, 10000. 

However, simply having a strong waiting-time depen- 
dence in a two-time correlation function does not imply 
glassy behaviour. Indeed, simple coarsening systems ex- 
hibit aging fll3|| To distinguish this type of dynamics, 
type I aging ]14[| , from glassy dynamics, type II aging, 
we have considered the overlap distribution function 



Q(t w + t, t w + t) 



{t + t w )sf\t + t w ). (5) 



This is measured by first relaxing the system from a disor- 
dered state for a time t w . Then two copies of the system 
are made, {S^} and and each is evolved for a 

further time t with independent thermal noise. 

In Fig. U we show a plot of Q(t w + t,t w + t) vs. 
C(t,t w ) for three different waiting times. In each case 
Q(t w + t,t w + t) — > as C(t,t w ) — * 0, indicating that 
the two copies of the system continue to move apart 
from each other irrespective of when the copies were 
made. This behaviour is indicative of type II, glassy ag- 
ing Note that this behaviour is not an automatic 
consequence of the existence of infinitely many ground 
states. For example, we expect that for a triangular lat- 
tice anti-ferromagnet, Q(t w + t,t w + t) would tend to 
a constant as C(t,t w ) tends to zero, despite the non- 
zero ground-state entropy per spin. Finally, we have 
observed that C(t, t w ) and Q(t w + t, t w + t) are related 
by C(2t,t w ) = Q(t w + t,t w + t) over a wide range of 
timescales. This relation is known to hold in equilibrium, 
where both functions are independent of t w |14|. It has 
also been verified in the aging regime of some trap mod- 



els as well as in the spin-glass version of the plaquette 
model 15 1. 

These results provide evidence that the free-energy 
landscape of the ferromagnetic plaquette model is in- 
deed complex and has many mutually inaccessible min- 
ima. This property is similar to what is seen in structural 
glasses and in mean field spin-glass models. 
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FIG. 5. Plot of Q(t w + t,t w + t) versus C(t,t w ). Three 
weighting times are shown: t w = 100(o), 1000(A), 5000(o). 



VI. CONCLUSIONS 



We have shown that the ferromagnetic plaquette model 
has many characteristics in common with glass forming 
liquids. We provided an explanation for the long-lived 
nature of the supercooled phase in terms of standard nu- 
cleation theory. By measuring the two-time correlation 
functions, we have observed two distinct kinds of behav- 
ior: (i) in the supercooled phase, the system appears to 
be stationary and the timescale of the energy autocorre- 
lation functions diverges at a temperature identified in 
earlier work as a glass transition temperature T g , and 
(ii) below T g the spin-spin autocorrelation functions ex- 
hibit aging of a form characteristic of glasses. Our find- 
ings, along with the work reported in ||,^]], indicate that 
a deeper understanding of this model would indeed be 
worthwhile. 
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